Reconstruction of the genome-scale metabolic network model of Sinorhizobium fredii CCBAU45436 for free-living and symbiotic states.

Sinorhizobium fredii CCBAU45436 is an excellent rhizobium that plays an important role in agricultural production. However, there still needs more comprehensive understanding of the metabolic system of S. fredii CCBAU45436, which hinders its application in agriculture. Therefore, based on the first-generation metabolic model iCC541 we developed a new genome-scale metabolic model iAQY970, which contains 970 genes, 1,052 reactions, 942 metabolites and is scored 89% in the MEMOTE test. Cell growth phenotype predicted by iAQY970 is 81.7% consistent with the experimental data. The results of mapping the proteome data under free-living and symbiosis conditions to the model showed that the biomass production rate in the logarithmic phase was faster than that in the stable phase, and the nitrogen fixation efficiency of rhizobia parasitized in cultivated soybean was higher than that in wild-type soybean, which was consistent with the actual situation. In the symbiotic condition, there are 184 genes that would affect growth, of which 94 are essential; In the free-living condition, there are 143 genes that influence growth, of which 78 are essential. Among them, 86 of the 94 essential genes in the symbiotic condition were consistent with the prediction of iCC541, and 44 essential genes were confirmed by literature information; meanwhile, 30 genes were identified by DEG and 33 genes were identified by Geptop. In addition, we extracted four key nitrogen fixation modules from the model and predicted that sulfite reductase (EC 1.8.7.1) and nitrogenase (EC 1.18.6.1) as the target enzymes to enhance nitrogen fixation by MOMA, which provided a potential focus for strain optimization. Through the comprehensive metabolic model, we can better understand the metabolic capabilities of S. fredii CCBAU45436 and make full use of it in the future.

Chromosome structure modeling tools and their evaluation in bacteria.

The three-dimensional (3D) structure of bacterial chromosomes is crucial for understanding chromosome function. With the growing availability of high-throughput chromosome conformation capture (3C/Hi-C) data, the 3D structure reconstruction algorithms have become powerful tools to study bacterial chromosome structure and function. It is highly desired to have a recommendation on the chromosome structure reconstruction tools to facilitate the prokaryotic 3D genomics. In this work, we review existing chromosome 3D structure reconstruction algorithms and classify them based on their underlying computational models into two categories: constraint-based modeling and thermodynamics-based modeling. We briefly compare these algorithms utilizing 3C/Hi-C datasets and fluorescence microscopy data obtained from Escherichia coli and Caulobacter crescentus, as well as simulated datasets. We discuss current challenges in the 3D reconstruction algorithms for bacterial chromosomes, primarily focusing on software usability. Finally, we briefly prospect future research directions for bacterial chromosome structure reconstruction algorithms.

EVRC: reconstruction of chromosome 3D structure models using error-vector resultant algorithm with clustering coefficient.

MOTIVATION: Reconstruction of 3D structure models is of great importance for the study of chromosome function. Software tools for this task are highly needed. RESULTS: We present a novel reconstruction algorithm, called EVRC, which utilizes co-clustering coefficients and error-vector resultant for chromosome 3D structure reconstruction. As an update of our previous EVR algorithm, EVRC now can deal with both single and multiple chromosomes in structure modeling. To evaluate the effectiveness and accuracy of the EVRC algorithm, we applied it to simulation datasets and real Hi-C datasets. The results show that the reconstructed structures have high similarity to the original/real structures, indicating the effectiveness and robustness of the EVRC algorithm. Furthermore, we applied the algorithm to the 3D conformation reconstruction of the wild-type and mutant Arabidopsis thaliana chromosomes and demonstrated the differences in structural characteristics between different chromosomes. We also accurately showed the conformational change in the centromere region of the mutant compared with the wild-type of Arabidopsis chromosome 1. Our EVRC algorithm is a valuable software tool for the field of chromatin structure reconstruction, and holds great promise for advancing our understanding on the chromosome functions. AVAILABILITY AND IMPLEMENTATION: The software is available at https://github.com/mbglab/EVRC.

Independent Component Analysis Reveals the Transcriptional Regulatory Modules in Bradyrhizobium diazoefficiens USDA110.

The dynamic adaptation of bacteria to environmental changes is achieved through the coordinated expression of many genes, which constitutes a transcriptional regulatory network (TRN). Bradyrhizobium diazoefficiens USDA110 is an important model strain for the study of symbiotic nitrogen fixation (SNF), and its SNF ability largely depends on the TRN. In this study, independent component analysis was applied to 226 high-quality gene expression profiles of B. diazoefficiens USDA110 microarray datasets, from which 64 iModulons were identified. Using these iModulons and their condition-specific activity levels, we (1) provided new insights into the connection between the FixLJ-FixK2-FixK1 regulatory cascade and quorum sensing, (2) discovered the independence of the FixLJ-FixK2-FixK1 and NifA/RpoN regulatory cascades in response to oxygen, (3) identified the FixLJ-FixK2 cascade as a mediator connecting the FixK2-2 iModulon and the Phenylalanine iModulon, (4) described the differential activation of iModulons in B. diazoefficiens USDA110 under different environmental conditions, and (5) proposed a notion of active-TRN based on the changes in iModulon activity to better illustrate the relationship between gene regulation and environmental condition. In sum, this research offered an iModulon-based TRN for B. diazoefficiens USDA110, which formed a foundation for comprehensively understanding the intricate transcriptional regulation during SNF.

OptoCRISPRi-HD: Engineering a Bacterial Green-Light-Activated CRISPRi System with a High Dynamic Range.

The ability to modulate gene expression is crucial for studying gene function and programming cell behaviors. Combining the reliability of CRISPRi and the precision of optogenetics, the optoCRISPRi technique is emerging as an advanced tool for live-cell gene regulation. Since previous versions of optoCRISPRi often exhibit no more than a 10-fold dynamic range due to the leakage activity, they are not suitable for targets that are sensitive to such leakage or critical for cell growth. Here, we describe a green-light-activated CRISPRi system with a high dynamic range (40 fold) and the flexibility of changing targets in Escherichia coli. Our optoCRISPRi-HD system can efficiently repress essential genes, nonessential genes, or inhibit the initiation of DNA replication. Providing a regulative system with high resolution over space-time and extensive targets, our study would facilitate further research involving complex gene networks, metabolic flux redirection, or bioprinting.

The Spatial Organization of Bacterial Transcriptional Regulatory Networks.

The transcriptional regulatory network (TRN) is the central pivot of a prokaryotic organism to receive, process and respond to internal and external environmental information. However, little is known about its spatial organization so far. In recent years, chromatin interaction data of bacteria such as Escherichia coli and Bacillus subtilis have been published, making it possible to study the spatial organization of bacterial transcriptional regulatory networks. By combining TRNs and chromatin interaction data of E. coli and B. subtilis, we explored the spatial organization characteristics of bacterial TRNs in many aspects such as regulation directions (positive and negative), central nodes (hubs, bottlenecks), hierarchical levels (top, middle, bottom) and network motifs (feed-forward loops and single input modules) of the TRNs and found that the bacterial TRNs have a variety of stable spatial organization features under different physiological conditions that may be closely related with biological functions. Our findings provided new insights into the connection between transcriptional regulation and the spatial organization of chromosome in bacteria and might serve as a factual foundation for trying spatial-distance-based gene circuit design in synthetic biology.

Defective chromatin architectures in embryonic stem cells derived from somatic cell nuclear transfer impair their differentiation potentials.

Nuclear transfer embryonic stem cells (ntESCs) hold enormous promise for individual-specific regenerative medicine. However, the chromatin states of ntESCs remain poorly characterized. In this study, we employed ATAC-seq and Hi-C techniques to explore the chromatin accessibility and three-dimensional (3D) genome organization of ntESCs. The results show that the chromatin accessibility and genome structures of somatic cells are re-arranged to ESC-like states overall in ntESCs, including compartments, topologically associating domains (TADs) and chromatin loops. However, compared to fertilized ESCs (fESCs), ntESCs show some abnormal openness and structures that have not been reprogrammed completely, which impair the differentiation potential of ntESCs. The histone modification H3K9me3 may be involved in abnormal structures in ntESCs, including incorrect compartment switches and incomplete TAD rebuilding. Moreover, ntESCs and iPSCs show high similarity in 3D genome structures, while a few differences are detected due to different somatic cell origins and reprogramming mechanisms. Through systematic analyses, our study provides a global view of chromatin accessibility and 3D genome organization in ntESCs, which can further facilitate the understanding of the similarities and differences between ntESCs and fESCs.

Computationally Reconstructed Interactome of Bradyrhizobium diazoefficiens USDA110 Reveals Novel Functional Modules and Protein Hubs for Symbiotic Nitrogen Fixation.

Symbiotic nitrogen fixation is an important part of the nitrogen biogeochemical cycles and the main nitrogen source of the biosphere. As a classical model system for symbiotic nitrogen fixation, rhizobium-legume systems have been studied elaborately for decades. Details about the molecular mechanisms of the communication and coordination between rhizobia and host plants is becoming clearer. For more systematic insights, there is an increasing demand for new studies integrating multiomics information. Here, we present a comprehensive computational framework integrating the reconstructed protein interactome of B. diazoefficiens USDA110 with its transcriptome and proteome data to study the complex protein-protein interaction (PPI) network involved in the symbiosis system. We reconstructed the interactome of B. diazoefficiens USDA110 by computational approaches. Based on the comparison of interactomes between B. diazoefficiens USDA110 and other rhizobia, we inferred that the slow growth of B. diazoefficiens USDA110 may be due to the requirement of more protein modifications, and we further identified 36 conserved functional PPI modules. Integrated with transcriptome and proteome data, interactomes representing free-living cell and symbiotic nitrogen-fixing (SNF) bacteroid were obtained. Based on the SNF interactome, a core-sub-PPI-network for symbiotic nitrogen fixation was determined and nine novel functional modules and eleven key protein hubs playing key roles in symbiosis were identified. The reconstructed interactome of B. diazoefficiens USDA110 may serve as a valuable reference for studying the mechanism underlying the SNF system of rhizobia and legumes.

EVR: reconstruction of bacterial chromosome 3D structure models using error-vector resultant algorithm.

BACKGROUND: More and more 3C/Hi-C experiments on prokaryotes have been published. However, most of the published modeling tools for chromosome 3D structures are targeting at eukaryotes. How to transform prokaryotic experimental chromosome interaction data into spatial structure models is an important task and in great need. RESULTS: We have developed a new reconstruction program for bacterial chromosome 3D structure models called EVR that exploits a simple Error-Vector Resultant (EVR) algorithm. This software tool is particularly optimized for the closed-loop structural features of prokaryotic chromosomes. The parallel implementation of the program can utilize the computing power of both multi-core CPUs and GPUs. CONCLUSIONS: EVR can be used to reconstruct the bacterial 3D chromosome structure based on the contact frequency matrix derived from 3C/Hi-C experimental data quickly and precisely.

Phosphates as Energy Sources to Expand Metabolic Networks.

Phosphates are essential for modern metabolisms. A recent study reported a phosphate-free metabolic network and suggested that thioesters, rather than phosphates, could alleviate thermodynamic bottlenecks of network expansion. As a result, it was considered that a phosphorus-independent metabolism could exist before the phosphate-based genetic coding system. To explore the origin of phosphorus-dependent metabolism, the present study constructs a protometabolic network that contains phosphates prebiotically available using computational systems biology approaches. It is found that some primitive phosphorylated intermediates could greatly alleviate thermodynamic bottlenecks of network expansion. Moreover, the phosphorus-dependent metabolic network exhibits several ancient features. Taken together, it is concluded that phosphates played a role as important as that of thioesters during the origin and evolution of metabolism. Both phosphorus and sulfur are speculated to be critical to the origin of life.

Spatial organization of the transcriptional regulatory network of Saccharomyces cerevisiae.

A transcriptional regulatory network (TRN) is a complex network composed of all of the regulatory interactions between transcription factors and the corresponding target genes. Recently, three-dimensional (3D) genomic studies have shown that the 3D structure of the genome may influence the regulation of gene transcription, which provides us with a novel perspective. In the present study, we constructed the TRN of the budding yeast Saccharomyces cerevisiae and placed it in the context of a 3D genome model. We analyzed the spatial organization of the yeast TRN on four levels: global features, central nodes, hierarchical structure and network motifs. The results obtained suggest that the TRN of S. cerevisiae presents an optimized structure in space to adapt to functional requirements.

Genomic analyses of primitive, wild and cultivated citrus provide insights into asexual reproduction.

The emergence of apomixis-the transition from sexual to asexual reproduction-is a prominent feature of modern citrus. Here we de novo sequenced and comprehensively studied the genomes of four representative citrus species. Additionally, we sequenced 100 accessions of primitive, wild and cultivated citrus. Comparative population analysis suggested that genomic regions harboring energy- and reproduction-associated genes are probably under selection in cultivated citrus. We also narrowed the genetic locus responsible for citrus polyembryony, a form of apomixis, to an 80-kb region containing 11 candidate genes. One of these, CitRWP, is expressed at higher levels in ovules of polyembryonic cultivars. We found a miniature inverted-repeat transposable element insertion in the promoter region of CitRWP that cosegregated with polyembryony. This study provides new insights into citrus apomixis and constitutes a promising resource for the mining of agriculturally important genes.

Construction and simulation of the Bradyrhizobium diazoefficiens USDA110 metabolic network: a comparison between free-living and symbiotic states.

Bradyrhizobium diazoefficiens is a rhizobium able to convert atmospheric nitrogen into ammonium by establishing mutualistic symbiosis with soybean. It has been recognized as an important parent strain for microbial agents and is widely applied in agricultural and environmental fields. In order to study the metabolic properties of symbiotic nitrogen fixation and the differences between a free-living cell and a symbiotic bacteroid, a genome-scale metabolic network of B. diazoefficiens USDA110 was constructed and analyzed. The metabolic network, iYY1101, contains 1031 reactions, 661 metabolites, and 1101 genes in total. Metabolic models reflecting free-living and symbiotic states were determined by defining the corresponding objective functions and substrate input sets, and were further constrained by high-throughput transcriptomic and proteomic data. Constraint-based flux analysis was used to compare the metabolic capacities and the effects on the metabolic targets of genes and reactions between the two physiological states. The results showed that a free-living rhizobium possesses a steady state flux distribution for sustaining a complex supply of biomass precursors while a symbiotic bacteroid maintains a relatively condensed one adapted to nitrogen-fixation. Our metabolic models may serve as a promising platform for better understanding the symbiotic nitrogen fixation of this species.

A novel proposal of a simplified bacterial gene set and the neo-construction of a general minimized metabolic network.

A minimal gene set (MGS) is critical for the assembly of a minimal artificial cell. We have developed a proposal of simplifying bacterial gene set to approximate a bacterial MGS by the following procedure. First, we base our simplified bacterial gene set (SBGS) on experimentally determined essential genes to ensure that the genes included in the SBGS are critical. Second, we introduced a half-retaining strategy to extract persistent essential genes to ensure stability. Third, we constructed a viable metabolic network to supplement SBGS. The proposed SBGS includes 327 genes and required 431 reactions. This report describes an SBGS that preserves both self-replication and self-maintenance systems. In the minimized metabolic network, we identified five novel hub metabolites and confirmed 20 known hubs. Highly essential genes were found to distribute the connecting metabolites into more reactions. Based on our SBGS, we expanded the pool of targets for designing broad-spectrum antibacterial drugs to reduce pathogen resistance. We also suggested a rough semi-de novo strategy to synthesize an artificial cell, with potential applications in industry.

Are protein hubs faster folders? Exploration based on Escherichia coli proteome.

Protein hubs in protein-protein interaction network are especially important due to their central roles in the entire network. Despite of their importance, the folding kinetics of hub proteins in comparison with non-hubs is still unknown. In this work, the folding rates for protein hubs and non-hubs were predicted and compared for the interactome of Escherichia coli K12, and the results showed that hub proteins fold faster than non-hub proteins. A possible explanation might be that protein hubs have more and fast-folding structural conformations than non-hubs, which leads to the notion of "hub of hubs" in the protein conformation space. It was found that the sequence and structure features relevant to protein folding rates are also different between hub and non-hub proteins. Moreover, the interacting proteins tend to have similar folding rates. These results gave insightful implications for understanding the interplay between the mechanisms of protein folding and interaction.

Abundance and Temperature Dependency of Protein-Protein Interaction Revealed by Interface Structure Analysis and Stability Evolution.

Protein complexes are major forms of protein-protein interactions and implement essential biological functions. The subunit interface in a protein complex is related to its thermostability. Though the roles of interface properties in thermal adaptation have been investigated for protein complexes, the relationship between the interface size and the expression level of the subunits remains unknown. In the present work, we studied this relationship and found a positive correlation in thermophiles rather than mesophiles. Moreover, we found that the protein interaction strength in complexes is not only temperature-dependent but also abundance-dependent. The underlying mechanism for the observed correlation was explored by simulating the evolution of protein interface stability, which highlights the avoidance of misinteraction. Our findings make more complete the picture of the mechanisms for protein complex thermal adaptation and provide new insights into the principles of protein-protein interactions.

Amino Acid Flux from Metabolic Network Benefits Protein Translation: the Role of Resource Availability.

Protein translation is a central step in gene expression and affected by many factors such as codon usage bias, mRNA folding energy and tRNA abundance. Despite intensive previous studies, how metabolic amino acid supply correlates with protein translation efficiency remains unknown. In this work, we estimated the amino acid flux from metabolic network for each protein in Escherichia coli and Saccharomyces cerevisiae by using Flux Balance Analysis. Integrated with the mRNA expression level, protein abundance and ribosome profiling data, we provided a detailed description of the role of amino acid supply in protein translation. Our results showed that amino acid supply positively correlates with translation efficiency and ribosome density. Moreover, with the rank-based regression model, we found that metabolic amino acid supply facilitates ribosome utilization. Based on the fact that the ribosome density change of well-amino-acid-supplied genes is smaller than poorly-amino-acid-supply genes under amino acid starvation, we reached the conclusion that amino acid supply may buffer ribosome density change against amino acid starvation and benefit maintaining a relatively stable translation environment. Our work provided new insights into the connection between metabolic amino acid supply and protein translation process by revealing a new regulation strategy that is dependent on resource availability.

Spatial features for Escherichia coli genome organization.

BACKGROUND: In bacterial genomes, the compactly encoded genes and operons are well organized, with genes in the same biological pathway or operons in the same regulon close to each other on the genome sequence. In addition, the linearly close genes have a higher probability of co-expression and their protein products tend to form protein-protein interactions. However, the organization features of bacterial genomes in a three-dimensional space remain elusive. The DNA interaction data of Escherichia coli, measured by the genome conformation capture (GCC) technique, have recently become available, which allowed us to investigate the spatial features of bacterial genome organization. RESULTS: By renormalizing the GCC data, we compared the interaction frequency of operon pairs in the same regulon with that of random operon pairs. The results showed that arrangements of operons in the E. coli genome tend to minimize the spatial distance between operons in the same regulon. A similar global organization feature exists for genes in biological pathways of E. coli. In addition, the genes close to each other spatially (even if they are far from each other on the genome sequence) tend to be co-expressed and form protein-protein interactions. These results provided new insights into the organization principles of bacterial genomes and support the notion of transcription factory. CONCLUSIONS: This study revealed the organization features of Escherichia coli genomic functional units in the 3D space and furthered our understanding of the link between the three-dimensional structure of chromosomes and biological function.

Sequence determinants of prokaryotic gene expression level under heat stress.

Prokaryotic gene expression is environment-dependent and temperature plays an important role in shaping the gene expression profile. Revealing the regulation mechanisms of gene expression pertaining to temperature has attracted tremendous efforts in recent years particularly owning to the yielding of transcriptome and proteome data by high-throughput techniques. However, most of the previous works concentrated on the characterization of the gene expression profile of individual organism and little effort has been made to disclose the commonality among organisms, especially for the gene sequence features. In this report, we collected the transcriptome and proteome data measured under heat stress condition from recently published literature and studied the sequence determinants for the expression level of heat-responsive genes on multiple layers. Our results showed that there indeed exist commonness and consistent patterns of the sequence features among organisms for the differentially expressed genes under heat stress condition. Some features are attributed to the requirement of thermostability while some are dominated by gene function. The revealed sequence determinants of bacterial gene expression level under heat stress complement the knowledge about the regulation factors of prokaryotic gene expression responding to the change of environmental conditions. Furthermore, comparisons to thermophilic adaption have been performed to reveal the similarity and dissimilarity of the sequence determinants for the response to heat stress and for the adaption to high habitat temperature, which elucidates the complex landscape of gene expression related to the same physical factor of temperature.